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Abstract 

A high reflectivity and polarization selective high contrast grating 
mirror has been designed with the use of an automated optimization 
algorithm. Through a precise study of the tolerance of the different 
lengths of the structure, the robustness with respect to the fabrication 
errors has been enhanced to high tolerance values between 5 % and 
210 %. This adjustment of the dimensions of the structure leads to a 
250 nm large bandwidth mirror well adapted for a VCSEL application 
at A = 2.65 |xm and can easily be scaled for other wavelengths. 
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1 Introduction 



The Tunable Diode Laser Absorption Spectroscopy (TDLAS) is a sensitive 
and fast method for gas sensing used for instance to measure concentra- 
tions of polluting gas such as CH4 and CO. This technique requires the use 
of a tunable, stable and single mode operating source emitting in the mid 
infrared wavelength range beyond 2 [im where strong absorption lines are 
present for these gas [H |2] . VCSELs are highly suitable for this kind of ap- 
plication, however, the development of continuous wave devices operating at 
room temperature at this wavelength range is still an important challenge. 

Due to the small thickness of the gain region, VCSEL structures require 
high quality cavities with highly reflective Bragg mirrors (R > 99 %). In 
the mid-infrared wavelength range, even though the best performances were 
obtained with AlGalnAsSb material system, the relatively low index con- 
trast (An ~ 0.5) between DBR layers make the mirrors become as thick as 
11 |j.m impairing the electro-thermal-optical properties of the structure [3]. 
Currently, laser emission has been shown from an all-epitaxial monolithic 
micro cavity near 2.3 (xm in a continuous wave mode and in quasi- CW (5 %, 
1 (is) up to 2.63 [im [4J. Another hybrid structure made of a dielectric top 
mirror and a buried tunnel junction operates in CW emission at room tem- 
perature in the 2.4-2.6 |am wavelength range 0. 

One solution to increase the emission wavelength is the use of a high 
contrast grating mirror as top cavity mirror [6j . Such gratings, combined with 
a low index sub-layer, can exhibit high reflectivity of more than 99.9 % for 
bandwidths larger than 100 nm [7]. Moreover, due to their one dimensional 
symmetry, a high polarization selectivity can be performed by these mirrors 
and with a total thickness of less than 2 ^m [S], the stability and quality of 
the VCSEL emitted beam should be improved. 

However, contrary to Bragg mirrors, the explanation of the optical re- 
sponse of high contrast grating (HCG) is not straightforward and the design 
adjustment made in order to achieve the required properties for a VCSEL 
application is still complex. Even though several formalism have been de- 
fined, based for instance on the destructive interference of modes P], the 
use of an optimization algorithm for the design combined with a numerical 
computation of the reflectivity keeps the advantage of versatility. Indeed, 
such a method easily allows the user to aim for specific properties of the 
HCG. It has been used for instance to design wide band and high diffraction 
efficiency grating [TU], large bandwidth grating mirror for both TE and TM 
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polarization [TT] or, like in this work, polarization selective large bandwidth 
mirror with technological constraints on the design dimensions [5]. 

In this work, we present a Si/Si02 HCG design optimized for a VCSEL 
application at A = 2.65 [J.m. Then, through a precise study of the computed 
tolerances of the structure dimensions, the robustness is enhanced with re- 
spect to the fabrication errors. 

2 Design of the HCG structure 

The high contrast grating structure studied in this work is based on Si/Si02 
materials which properties and fabrication process are well known. The grat- 
ing is made of silicium (n = 3.435) on top of a low index layer of Si02 
(n = 1.509) used to achieved high reflectivity [7] and allowing the use of a 
selective etching method for the fabrication of the grating. In order to in- 
crease the total reflectivity of the reflector and broaden the stopband, two 
quarter wavelength layers of Si/Si02 are combined with the grating [S]. For 
the simulation of the design, the substrate is chosen as the VCSEL cavity 
material with an optical index of n = 3.521. 

In order to be well suitable for a VCSEL application, HCGs have to ex- 
hibit optical properties as we have defined with a 99.9 % transverse magnetic 
(TM) reflectivity for the largest possible bandwidth. Moreover, to ensure a 
polarization stability of the emitted beam, the reflectivity of the transverse 
electric (TE) mode has been chosen to be kept lower than 90 % for the whole 
bandwidth. 

A 99.5 % TM reflectivity should be enough to achieve laser emission but 
a 0.4 % security margin has been chosen to account for possible experimen- 
tal growth imperfections and losses due to absorption. Indeed, the presence 
of OH radicals as impurities in Si02 results in an absorption band in the 
2.6-2.9 [im range [12]. In this reference, the absorption of silica has been 
measured with a 10 dB/cm value at 2.65 |am, i.e. the refractive index be- 
comes nsi02 = 1-509 + 2.1e~^j, which results in a 0.1 % fall of the HCG 
reflectivity while a 0.4 % decrease of reflectivity has been observed for a sim- 
ulation with a 20 dB/cm attenuation coefficient ( nsi02 = 1-509 + 1.7e^^j ) 
at 2.675 \xm. This latter value appears therefore as the maximum absorption 
value allowed for the studied structure. However, the OH impurities con- 
centration should strongly depends on the fabrication process and the exact 
value of the absorption must be determined in each case. For this theoretical 
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Figure 1: Scheme of the structure. The sihcium grating defined by the empty 
length of the grooves Lg, the filled length Lf and the grating thickness Tg is 
on top of a Si02 layer of thickness T^. These four parameters are adjusted 
to meet the characteristics of a VCSEL refiector. 

study, a pure Si02 without any OH absorption was considered. Nevertheless, 
if the silica absorption is proven to be too high, this material can be replaced 
by another dielectric such as Si3N4. 

The evaluation of the refiectivity of the mirror is made thanks to a rigor- 
ous coupled wave analysis (RCWA) [131 E] which numerically finds an exact 
solution of Maxwell's equations for the electromagnetic diffraction problem 
of an infinite grating structure. 

Refiection spectra for TE and TM modes are thus computed and through 
a well defined quality factor [8], the optical performance of the refiector is 
numerically evaluated. This performance is then maximized by a genetic 
optimization algorithm which searches for a global maximum of the quality 
factor by adjusting the design. The use of a global algorithm is manda- 
tory in this problem since the quality factor function exhibits numerous local 
maxima. 

In this work, the performances of the grating are optimized by adjusting 
the empty and filled lengths (Lg and Lf), the grating thickness Tg and the 
sublayer thickness as shown on Figure [TJ 

The use of an automated optimization presents the advantage of easily 
imposing technological constraints on the dimensions of the structure. These 
constraints have been set on the lengths Lg and Lf limited to a minimum 
value of 500 nm which should ease the photolithographic process. Besides, 
the etching process is also a critical step of the grating fabrication and, to 
achieve the theoretical grating profile presented on Figure [H the shape factor 
SF = Le/Tg of the grooves is kept at a minimum value of 0.9 since squared 
patterns are easier to etch than deep ones. 
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Figure 2: Reflection spectra for tlie TM mode (blue) and TE mode (dashed 
red) of the structure automatically optimized by a genetic-based algorithm. 
This optimum design exhibits a 307 nm large bandwidth with a 99.9 % high 
TM reflectivity. 

The optimization for A = 2.65 |am of the structure shown on Figure 
[U results in a 307 nm large bandwidth mirror exhibiting a high polariza- 
tion selectivity with Rte < 90 % (Fig. [2]). The optimum lengths of the 
design are given in Table [1] and satisfy the technological constraints with 
Lg = 799 nm, Lf = 541 nm and a shape factor SF = 0.9. The grating 
thickness Tg = 886 nm and the sublayer thickness = 321 nm combined 
with the two quarter wavelength layers result in a 1.84 [j.m thick reflector. 
Even though the performances of this HCG are well adapted for a VCSEL 
application, the tolerances of the dimensions of this grating are critical. For 
instance, the 884 nm minimum grating thickness leads to a maximum error 
allowed as small as 2 nm. 

3 Optimization of the robustness 

The evaluation of the tolerances is performed by computing the variation 
range of one parameter for which the HCG keeps a 99.9 % TM reflectivity 
together with a Rte < 90 % at Aq = 2.65 \im. This computation is made by 
varying one parameter at a time, for instance the grating thickness Tg, while 
keeping the other ones (T^, and Lf) constant at their optimal value. 

The computation of the tolerances of the design found by the optimiza- 
tion algorithm exhibits large variation ranges of ALg = 14 %, ALj = 25 %, 
ATg = 5 % and AT^ = 68 %. However, the optimum point is not centred in 
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Table 1: Tolerances of the optimum design found by the genetic-based algo- 
rithm. 







Optimum 


Minimum 


Maximum 


Lf 
T, 
Tl 




799 nm 
541 nm 

886 nm 
321 nm 


689 nm 
495 nm 

884 nm 
196 nm 


803 nm 
629 nm 

932 nm 
414 nm 


A = 
FF 


Le + Lf 

= Lf/A 


1340 nm 
40.37 % 


1260 nm 
35.97 % 


1345 nm 
45.67 % 


de = 

/3e = 


--L,+Tg 

-- Le/ae 


1685 nm 
52.58 % 


1652 nm 
50.80 % 


1689 nm 
54.89 % 


«/ = 
Pf = 


-Lf + Tg 
-- Lf/af 


1427 nm 
62.09 % 


1388 nm 
59.13 % 


1430 nm 
64.74 % 



these variation ranges and the real tolerance values are much more limited 
with for instance a maximum increase of 4 nm on the empty length Lg. More- 
over, these variation ranges are computed separately and give no information 
of the simultaneous error allowed on several parameters. 

A solution to access this information is to compute the variation range 
of combinations of two dimensions. In the following, only {Lg, Lf), {Lg, Tg) 
and (L/, Tg) combinations are evaluated since the sublayer thickness Tl is 
the most tolerant and centred parameter. Besides, Tl does not depend on 
the etching process and should be fabricated with a better accuracy than the 
grating parameters. 

For the first couple (Le, I//), the variation ranges of two significant com- 
binations of Le and Lf are evaluated. The first one is the grating period 
given hy K = Le + Lf. This parameter is varied around its optimum vahie 
of A = 1340 nm while keeping the fill factor defined as FF = Le/A at its 
optimum value of 40.37 %. The second combination of Lg and Lf is the fill 
factor FF which variation range is evaluated with a constant grating period 
of A =1340 nm. Large tolerance values are exhibited by these parameters 
with 85 nm for AA and 9.7 % for AFF. The minimum and maximum of the 
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Lf (nm) 




o Optimum (799, 541) 
★ Centre (753, 560) 
L^i„=(689,541) 
La^ax=(803, 541) 
Lf,„in= (799, 495) 



Lp (nm) 



= (799, 629) 
= (858, 481) 
= (728, 611) 
(751, 508) 
: (801, 543) 



Figure 3: Tolerance map of Lg and Lf. The variation ranges of the grating 
period A = Le + Lf and the fill factor FF = Lf/A define a polygon (grey) 
of allowed (Lg, Lf) couples for the design. The centre of the incircle {-k) 
enhances the robustness of the optimum point (o). 

variation ranges of the four parameters Lg, Lf, A and FF, summarized in 
Table [H creating eight (Lg, Lf) couples which can be plotted in a. Lf versus 
Le graph (Fig. [3]). In the (Lg, Lf) plan, these points define a polygon, in 
grey on Figure [3l which area represents the tolerance area of (Lg, Lf) couples 
allowed for the design. 

The representation of the tolerances in a (Lg, Lf) plan also shows very 
well the position of the optimum point within the tolerance area. In this case, 
the optimum point (o on Figure [3]) is located at an edge of the tolerance area 
(Lg = 799 nm, Lf = 541 nm). In order to increase the robustness regarding 
to the errors of fabrication which could be made on Lg and Lf, the design 
should be centred within the tolerance area. To do so, the centre of the 
largest incircle of the polygon representing the tolerance area is computed. 
This point is thus the farthest point from any edge of the polygon. The 
centre (*) on Figure [3] corresponds to Lg = 753 nm and Lf = 560 nm and 
with a incircle radius of 47 nm, ensures a minimum tolerance of 94 nm on 
any combination of (Lg, Lf). 

The second dimension couple studied is the empty length and grating 
thickness (Lg, Tg). In this case, both values of the optimum point (799, 886) 
are closed to the tolerance limit (Lg max = 803 nm, Tg min = 884 nm). This 
can easily be seen on Figure H] where the optimum is localized at the edge of 
the tolerance area. In this map, the first combination of (Lg, Tg) is defined 
by Qfg = Lg + Tg and the second one by /3g = L^/ae. The polygon defined 
by the height extrema of the variation ranges of Lg, Tg, and /3g (Tabled]) 
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Tg (nm) 




o Optimum (799, 886) 
■k Centre (772, 898) 



(689, 886 



La„ax= (803, 886 
Vm= (884,8 
Tg™x= (932,886 
Pemin= (829,855 

Pemax= (760, 924 

aemin= (783, 868 

aerr,ax= (800, 888 



Figure 4: Tolerance map of Le and Tg. The variation ranges of Oe = Le + Tg 
and /3e = Le/ae define a polygon (grey) of allowed (Le, Tg) couples for the 
design. The centre of the incircle {-k) enhances the robustness of the optimum 
point (o). 



oOptimum (541, 886) 
★ Centre (519, 894) 
Lfmin= (495, 886) 
L,rr,ax= (629, 886) 
Tip,m= (541, 884) 
Tip,ax= (541, 932) 
Pf^i„= (583, 843) 
Pf™x= (503,923) 
a,min= (526, 861) 
If max = (542, 887) 



Figure 5: Tolerance map of Lf and Tg. The variation ranges of Of = Lf + Tg 
and (3f = Lf/af define a polygon (grey) of allowed {Lf, Tg) couples for the 
design. The centre of the incircle (*) enhances the robustness of the optimum 
point (o). 

exhibits an incircle with a radius of 28 nm centred at (772, 898) making a 
tolerance of a minimum value of 56 nm on the lengths (Lg, Tg) and any kind 
of their combinations. 

The last couple, formed by the filled length Lf and grating thickness 
Tg, is optimized by computing the tolerance values of af = Lf + Tg and 
/3f = Lf /af. As shown on Figure [5l the polygon exhibits two large areas 
joined by a very thin path where the optimum point is located at (541, 886). 
Once again, the optimum point is in a critical location closed to several 
edges of the tolerance area. By choosing the centre of the incircle of the 
polygon, located at (519, 894), variation ranges of the lengths Lf and Tg can 
be increased to a 44 nm value. 
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Table 2: Tolerances of the resulting design with optimized variation ranges. 





Optimum 


Minimum 


Maximum 






773 nm 


703 nm 


821 nm 


(15 %) 


Lf 


541 nm 


504 nm 


608 nm 


(19 %) 


T, 


894 nm 


870 nm 


924 nm 


(6 %) 


Tl 


321 nm 


nm 


674 nm 


(210 %) 




Wavelength (|jm) 



Figure 6: Reflection spectra for the TM mode (blue) and TE mode (dashed 
red) of the robust design with optimized tolerance values. This mirror ex- 
hibits a 250 nm large bandwidth for Rtm > 99.9 % and between 5 % and 
210 % of tolerance on its dimensions. 

4 Characteristics of the robust HCG 

As a result of this optimization of the tolerances, three different designs are 
obtained. Each one is optimized to have the best tolerances for the cou- 
ples (Le, Lf), (Le, Tg) and {Lf, Tg). These three points represent a triangle 
in a three dimensional space (Le, Lf, Tg). Thus, to flnd a unique design 
with optimized tolerances, the centre of the incircle of this triangle has been 
computed and result in a flnal optimized design with lengths Le = 773 nm, 
Lf = 541 nm and Tg = 894 nm. This structure exhibits a 250 nm large 
bandwidth as shown on Figure Ej which is 59 nm less than the optimum one 
found by the genetic-based algorithm. However, the computation of the vari- 
ation ranges of the design dimensions (Table [2]) exhibits a very robust design 
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with ALe = 15%, ALf = 19%, ATg = 6% and ATl = 210%. Each optimum 
dimension of the structure is well centred within the variation ranges which 
would ease the fabrication process. 



5 Conclusion 

In this work, a Si/Si02 high contrast grating mirror has been designed for a 
VCSEL application at 2.65 [im. This mirror exhibits a 250 nm large band- 
width for Rtm > 99.9 % together with a strong polarization selectivity by 
keeping Rte < 90 %. Moreover, with a total thickness of less than 2 |am, 
such reflector should improve the quality of the emitted laser beam of the 
VCSEL. 

From a technological point of view, the fabrication constraints are re- 
spected with a large pattern resolution (> 500 nm) and a grating shape 
factor closed to 0.9. Moreover, the robustness with respect to the fabrication 
errors has been enhanced and leads to tolerance values larger than 5 % on the 
structure dimensions. Such characteristics make this HCG well adapted for 
a VCSEL application and should limit the pitfalls during the manufacturing 
process. Moreover, as spotted in |7], HCG can be scaled with wavelength in 
the limit of the refractive index dispersion. 
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